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Abstract 

As a consequence of the LSND anomaly and other hints of an eV scale sterile neutrino from 
particle physics and cosmology, the neutrino sector of the standard model of particle physics has 
to be extended and the smallest extension is the (3+1) model, i.e., three active neutrinos plus one 
sterile one. In this work we study the neutrino mass matrix M u with three texture zeros in the (3+1) 
model, assuming all the neutrino states are of Majorana type. With the mass hierarchy between the 
active and sterile neutrinos as well as the smallness of the reactor neutrino mixing and active-sterile 
mixing, the analytical expressions can be greatly simplified. We systematically examine all the 120 
texture zeros, via both analytical and numerical analysis, and find that the 100 textures with zeros 
in the fourth column and row of M v can hardly be compatible with experiments; whereas, 19 of 
the other 20 are feasible. We can foresee that the active neutrino masses are required to be around 
0.01 eV and the next generation neutrinoless double beta decay experiments are promising to test 
these textures. Many other phenomenologically meaningful predictions are also obtained. All the 
CP conserving (3+1) textures with three zeros are excluded by experiments. 
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I. INTRODUCTION 



After observation of the smallest neutrino oscillation by Daya Bay [1[ for the first time 
at the CL. of 5cr, and immediately further confirmed by RENO |2|, dozens of solar, at- 
mospheric, reactor, and accelerator-based neutrino experiments have compellingly provided 
us with solid information on five physical parameters involving three active flavors of mas- 
sive neutrinos The three mixing angles quantify the difference between the flavor and 
mass eigenstates of the three generations of neutrinos and the amplitude of neutrino oscilla- 
tions; the two independent mass squared differences are directly connected to neutrino mass 
spectrum and crucial to the oscillation frequencies in neutrino experiments. Although our 
knowledge on neutrinos in the standard model (SM) of particle physics has achieved great 
progress, many basic properties of these magic particles are unclear and to be explored, e.g., 
their absolute masses, the mechanism for their tiny masses, and CP violation in the lepton 
sector. 

One puzzle on neutrinos is the long-standing LSND anomaly 0], which points to the exis- 
tence of one eV scale neutrino and cannot be accommodated in the three-neutrino paradigm. 
The anomalous v v — > v e in LSND motivated the MiniBooNE experiment whose latest low 
energy data support the LSND result, although the neutrino and antineutrino mode of the 
MiniBooNE experiment are incompatible to a large extent 0] . Furthermore, recalculation of 
the reactor v e flux indicates the anomalous v e disappearance 0], i.e., the so-called "reactor 
anomaly". The Gallium radioactive source experiments GALLEX and SAGE [8[ have also 
reported hints of v e disappearance, corresponding to the eV scale neutrino oscillation. In 
addition, cosmological data on cosmological microwave background, large scale structure, 
and the Hubble constant imply extra light degrees of freedom j^] , and light extra neutrinos 
beyond the three in the SM are convincingly the most natural candidate. All these hints of 
extra neutrinos, however, are restricted by, e.g., the ICARUS and v e — 12 C scattering [111 
data. 



Due to the stringent limit on the Z boson decay width [12j, the extra light neutrinos 
cannot directly participate in the weak interaction and is consequently sterile. The simplest 
extension of the lepton sector is the addition of one sterile flavor to the three active ones, i.e., 
the resultant (3+1) neutrino model. In this paper we concentrate on the simplest extension, 
although introducing more sterile states may reconcile comparatively better the data sets 



on sterile neutrinos and obtain a more preferable goodness-of-fit |13l . [14 . 

Neutrino mass matrix M v with texture zeros, i.e., vanishing elements, is an attractive 
framework to phenomenologically study neutrinos and the connection of neutrino physics 
to SM and new physics beyond SM. Texture zero structures imply the underlying flavor 
symmetry or mechanism [l5|, e.g., the Z n symmetry or Froggatt-Nielsen mechanism (l6| . 
Furthermore, with vanishing elements in M u , the physical parameters on neutrino mass, 
mixing and CP violation are always connected to each other. The framework of texture 
zeros have been intensively studied in the three-neutrino frame and it can likewise 

apply to the (3+1) model with an extra sterile neutrino 25j. In the (3+1) model, if all the 
four neutrino states are of Majorana type, the neutrino mass matrix M v has ten independent 
elements, 

/ m ee m e/i m er m es \ 



M„ 



Tflefj, ^fif^i TTlfiT TTlfis 
Tfl er Tfl^q- Tfl TT Tfi TS 

f^es TTlus T^Its 7^1 ss 
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compared to six in the three-neutrino scheme. Ghosh et al. 25| has recently investigated 
the (3+1) model with two zeros, analyzing systematically all the C 2 = 45 distinct textures. 
They found that the 15 textures with both the two zeros in the active sector (the first three 
rows and columns) of M v are compatible with current neutrino data, whereas the 30 with 
zero(s) in the sterile sector (the fourth row and column) are excluded by experiments mainly 
due to the mass hierarchy between the active and sterile neutrinos. 1 

In the presence of CP violation phases, one zero in M v generally imposes two correlation 
conditions on the parameters on neutrino masses, mixing angles and CP violating phases, 
i.e., one is related to the vanishing absolute value and another to the vanishing argument. 
In general, more zeros means more stringent constraints on neutrino parameters and conse- 
quently more accurate predictions. In this work we focuses on 4 x 4 neutrino mass matrices 
with three vanishing elements, assuming all the four neutrino states are of Majorana type. 
As a comparison, none of the 20 3 x 3 mass matrices with three zeros are compatible with 
current data 17J, |l8j. Adding one sterile neutrino, the number of possible textures with 
three zeros sum up to Cf = 120 and this work aims to investigate the feasibility of three 
zero textures in the (3+1) scheme and explore their phenomenological implications. 

Some elements of the lepton mixing matrix are more lengthy in the (3+1) model than in 
the three neutrino scheme [27] and resultantly most of the elements of M v become even more 
lengthy. However, due to the large mass hierarchy between the sterile and active neutrinos 
m4 3> 77ii,2,3 and the smallness of the angles governing reactor mixing [l], 0] and active- 



sterile mixing (ASM) [13|, |14|, l28H30j . the vanishing conditions m a p = (a, (3 = e, /i, r, s) 
can be unbelievably simplified and we arrive at clean expressions for the three ASM angles. 
Based on these expressions we predict the restricted ranges of the ASM angles and other 
phenomenological observables, e.g., the effective Majorana electron neutrino mass (m ee ) in 
neutrinoless double beta decay (0v(3f3). 

We also numerically analyse these textures and study their phenomenological implica- 
tions. We apply experimental constraints on the five active neutrino oscillation parameters, 
making use of the 3a ranges from Ref. 31], and examine which textures can accommodate 



the eV scale sterile neutrino mass and ~0.1 ASMs implied by short baseline (SBL) experi- 
ments. The sterile neutrino fits [13L 14, 28-301 are somehow different from each other due 



to variations of the experimental data and calculation procedure, and we assign the fit in 
Ref. ' 



29] as the criterion on the three parameters: 



LSND mass square difference Am^, the 



ASM angle 644 relating v e and u s , and #24 relating v ^ and v s . Specifically, as cosmological 
data prefer lighter neutrino states, we use the 3a parameter regions in the Am| x — sin 2 26 ee 
and Am^j — sin 2 28^ planes with ~ eV 2 Am^ in Fig. 1 of Ref. j29|, where the effective 
mixing angles are defined as 



sin 



29 r 



41 U. 



a l 



\u, 



a 1 



(2) 



with a = e, fi. SBL experiments are not sensitive to the neutrino of tauon flavor, then the 



constraint on the third ASM angle #34 is comparatively rather loose, 0° < 634 < 30° [26 



1 Actually their analysis is inadequate and not all the 30 structures with zero(s) in the sterile sector are 
definitely excluded. If m TS , as well as one matrix element in the active sector, vanishes, one resultant 
possibility is that the mixing angle #34 is close to zero. As the short baseline experiments are not sensitive 



to the tauon flavor, the vanishing of 634 is allowed [26| and it is yet immature to assert that such textures 



are excluded. See the following analytical and numerical discussions for details. 
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In our numerical calculation, all the randomly generated parameters constrained by exper- 
iments are evenly distributed in their 3a ranges. Other parameters, e.g., the CP violating 
phases £13, have not been determined at the 3a C.L. If they are randomly generated, we 
assume they distribute evenly in their logically allowed ranges, e.g., 5 i3 G [0, 2ir]. 

The paper is organized as follows: The 120 textures with three zeros in the (3+1) model 
are systematically classified and labeled in the Appendix. The formulas and experimental 
constraints on the neutrino parameters are briefly exposited in Sec. [TT1 Section HTT1 is devoted 
to analytical analysis of the textures and Sec. IIVI to the numerical discussions. We briefly 
comment on the textures without CP violation in Sec. |V]and end up with the conclusion in 
Sec.|VTl 



II. FORMULAS AND EXPERIMENTAL CONSTRAINTS 

Before demonstration of the analytical and numerical calculations, it is necessary to clarify 
our notations of the neutrino parameters and briefly explain relevant formulas, which is the 
aim of this section. We also sketch the current experimental constraints on the neutrino 
parameters. 



Notations and formulas 



Assuming the charged lepton mass matrix is diagonal and positive-definite, the 4x4 
Majorana neutrino mass matrix M v can be expressed as, in the flavor basis, 



(3) 



where M^ iag = diagjmi, m2, m^, m^} is the diagonal neutrino matrix in the mass basis, and 
V is the 4x4 PMNS lepton mixing matrix. In general, V can split into two parts, 



V = UP 



(4) 



U is the lepton mixing matrix for Dirac neutrinos, containing the solar mixing angle #12, 
atmospheric mixing angle 6*23, reactor mixing angle 9 13 , three ASM angles 6 ) 14 , #24, and #34, 



and three Dirac CP violating phases. If we parametrize U explicitly as [32 

U = (-R34-R24-R14) (-R23-R13) R12 , 



(5) 



the three Dirac phases are 613, 814, and 624 , where Rij/Rij is the rotation matrix in the ij 
flavor space, e.g., 
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(6a) 



(6b) 
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with Sij 



sin 6ij and 



matrix written as 



cos 9ij. The matrix P in Eq. (j4j) is a diagonal Majorana phase 
P = diag{l,e iQ/2 ,e i/3/2 ,e^ /2 }. (7) 



The factors 1/2's are introduced as a consequence of the mixing matrix V appearing twice 
in Eq. 

Noticing that M^ iag is sandwiched by the diagonal phase matrix P, we define 

(8a) 
(8b) 
(8c) 
(8d) 



Ai 


= m x , 


A 2 


= m 2 e ia 


A 3 


= m 3 e l/9 


A 4 


= m 4 e 17 



and the 4x4 neutrino mass matrix is slightly simplified to 

M v = U diag{Ai, A 2 , A 3 , A 4 } U T . 
When three independent entries of M v are set to zero, 



M^ ab 
M„ >cd 

M v ,ef 



m ab ■ 
m cd ■ 
m ef 



0, 
0, 

: 



where a, b, c, d, e, / = e, /i, r, s, or, explicitly, 

XiU al U bl + x 2 u a2 u b2 + x 3 u a3 u b3 + x 4 u a4 u M = , 

XiU cl U dl + X 2 U c2 U d2 + X 3 U c3 U d3 + X 4 U c4 U d4 = , 
XxUexUf! + X 2 U e2 U f2 + X 3 U e3 U f3 + X A U e4 U fA = 



(9) 

(10a) 
(10b) 
(10c) 



(11a) 
(lib) 
(11c) 



As we have only three complex equations, less than the number of neutrino flavors in the 
(3+1) model, we divide A 2 ,3,4 by Xi = mi and obtain the three following quantities: 



A21 
A31 
A41 



A2/A1 = — e l 
mi 



A3/A1 



mi 

A 4 /Ai = !^ e *r . 
mi 



(12a) 
(12b) 
(12c) 



With three zeros in the 4x4 matrix M u , based on Eqs. ffTTl) and f|T2|) . A 2 i i3 i j 4i can be 
expressed as functions of the matrix elements of U, 



A21 
A 3 i 
A41 



A2IJV/ A 23 41 j 

X 3 in/ X 234 i , 

A4IJV/ A2 3 41 , 



where 



A 



23 11 



U a 2U b2 U a3 U b3 U a4 U b4 
U c2 U d2 U c3 U d3 U c 4U d 4 

u e2 u f2 U e3 U f3 U e4 U f4 



(13a) 
(13b) 
(13c) 



(14a) 
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217V 



A 



317V 



A 



417V 



UalUbl U a3 Ub3 U a4 UbA 
UdUdl U c sUd3 UdUdA 

u el u fl U e3 U f3 U e4 U f4 

U a lUbl U a2 Ub2 U a JJb1 
UdUdl U c2 U d2 UctPdA 

UciU S i u e2 u f2 u ei u f4 

UalUbl U a2 Ub 2 U a3 Ub3 

UdUdi U c2 U d2 Uc 3 U d3 
UciU S i U e2 U f2 U e3 U f3 



(14b) 
(14c) 
(14d) 



In presence of the three Dirac CP violating phases S13, #14 and S 2 ^, the matrix U is generally 
complex, then, from the definitions of A2i,3i,4i, we can extract three independent neutrino 
mass ratios 



(15a) 
(15b) 
(15c) 



m 2 


A21 


= A 2 itv/ A 


mi 




m 3 _ 


A31 


= A 3 itv/ A 


mi 




7724 


A41 


= \^4in/ A 


mi 





and the three Majorana CP violating phases 

a = arg(A 2 i) 
(3 = arg(A 3i ) 
7 = arg(A 4i ) 



arg(A 2 i7v/A 2 34i) , 

arg(A 3 iAr/A234l) , 
arg(A 4 iAr/A234l) • 



(16a) 
(16b) 
(16c) 



Given the three vanishing conditions on M u , i.e., Eq. f llOp . the neutrino mass ratios and 
Majorana phases are related to the six lepton mixing angles and the three Dirac CP violating 
phases. To further relate these neutrino parameters to neutrino oscillation experiments and 
utilize the experimental constraints on the solar neutrino mass square difference Am^, the 
atmospheric neutrino mass square difference (Am^l as well as Am^, we define the two 
independent neutrino mass square difference ratios 

Ryl 
Ru2 



Am 2 2 i 




A21 


2 -l 


\Amli 




A31 


2 - 1| 



Amli 



I A 



21 



I A 



(17a) 
(17b) 



B. Experimental constraints on neutrino parameters 

For the convenience of the analytical and numerical calculations below, let us first clearly 
count the relevant neutrino parameters we are confront with. In presence of the extra sterile 
flavor, we have altogether 16 physical parameters 

• Six mixing angles between different flavors: three denote the mixing between active 
flavors, #12,23,13; while the other three are relevant to the ASMs, #14,24,34- 
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TABLE I: Experimental data on neutrino oscillation parameters 31] . 



Parameter 


Hierarchy Best fit 


3a 


range 


Ai?4 [l(T 5 eV 2 ] 


NH/IH 


7.62 


7.12 


- 8.20 


lAmfil [10~ 3 eV 2 ] 


NH 


2.55 


2.31 


- 2.74 




IH 


2.43 


2.21 


- 2.64 


#12 


NH/IH 


34.4° 


31.3° 


- 37.5° 


#23 


NH 


51.5° 


36.9° 


- 55.6° 




IH 


50.8° 


37.5° 


- 54.9° 




NH 


9.0° 


7.5° 


- 10.5° 




IH 


9.1° 


7.5° 


- 10.5° 



• Six CP violating phases: three are of Dirac type, #13,14, 24 j and the other three are of 
Majorana type, a, (3, and 7. 

• Four neutrino masses: m 1 , m 2 , m 3 , and m 4 , with the first three predominantly of 
active flavors, and the last one predominantly of sterile. 

Currently, the three active neutrino mixing angles and two mass squared differences Am^ 
and jAm^l can be extracted from neutrino oscillation data. Their corresponding central 
values and 3a ranges are listed in Table [B The Dirac CP violating phase #13 in the (3+1) 
model corresponds to the sole Dirac phase 5 in the three-neutrino scheme. Although it can 
be loosely constrained in global neutrino fits, we treat it as a free parameter due to the low 



sensitivity of current data to CP violation [3l|, [33 



We assume that the sterile neutrino is heavier than the three active ones, i.e., the LSND 
mass squared difference is positive Am^ > 0. In the opposite case Am^ < with three 
heavier active neutrinos, the sum of neutrino masses E m will be ~ 3 eV, which is much 
less favored than the former case by cosmological observations 0]. The sign of Am^ has 
not been pinned down and then two mass spectra are allowed: normal hierarchy (NH) with 
r«3 > mi, 2 and inverted hierarchy (IH) with 7713 < mi, 2- The lightest neutrino mass mo in 
the two spectrums are, respectively, rriQ = mi and mo = m^. 

As mentioned in the introduction, the 3a bound on Am^ and #14 04 used in this work 



is from Fig. 1 of Ref. [29] and the constraint on #34 is from Ref. [26], as presented in 
Table [III Beware that the upper and lower bounds of 6*14 and 6*24 are dependent of Am^ 
and the bounds given in Table [TT] are only the allowed maximum and minimum values. 
It is phenomenologically meaningful to investigate whether these three-zero textures are 
compatible with the eV sterile neutrino mass and small ASM. To this end, Am^ treated 
as input parameters, we calculate analytically and numerically predictions of the textures 
on the three ASM angles and compare them with experimental constraints. Consequently, 
if the 16 neutrino parameters are roughly classified into the "know" and the "unknown", 
we have six "known" parameters: the five active neutrino parameters (#12,23,135 Am^, and 
lAmJ-jJ) and Am^. The left 10 are otherwise treated as "unknown": the three ASM angles, 
the six CP violating phases, and the lightest neutrino mass m . 

In general, the three vanishing elements of M v induce six independent relations among the 
neutrino masses, mixing angles, and CP violating phases, i.e., Eq. f fl5|) and ffl~6]) . That is to 
say, the 10 "unknown" neutrino parameters are constrained by six correlation conditions and 
we are left with four degrees of freedom completely unconstrained. In other words, generally 
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TABLE II: Current constraints on sterile neutrino parameters 26|, |29(]. See text for the comment 
on the angles #14,24- 



Parameter Hierarchy 


3a 


range 


Am 2 41 [eV 2 ] NH/IH 


0.72 


- 2.53 


014 NH/IH 


6.1° 


- 15.3° 


6*24 NH/IH 


2.8° 


- 12.5° 


<9 34 NH/IH 


0° 


- 30° 



we cannot completely determine the 10 "unknown" parameters in the (3+1) model with 
three zeros. However, according to the analytical and numerical calculations below, we are 
able to obtain phenomenologically meaningful predictions on neutrino mass spectrum, ASMs 
and CP violating phases. Furthermore, we can also foresee the constrained ranges of the 
nonoscillation quantities: the sum of neutrino masses 

4 

S ro = ^ rxij , (18) 

i=l 

the effective electron neutrino mass in tritium beta decay 



4 

1=1 

and the effective Majorana electron neutrino mass governing Ovfifi: 

(20) 

One may notice that if all the three Dirac CP violating phases vanish, according to the 
relations in Eq. ( fl"6j) . the three Majorana phases are also zero, and the lepton mixing matrix 
U is real. 2 In this simplified case, we are left with only three correlation conditions of form 

= Ail = AjiAr/A2341 , (21) 

mi 

with i — 2, 3, 4, and the ten "unknown" parameters shrink to four, i.e., the three ASM angles 
and the lightest neutrino mass rriQ. Clearly, the correlation conditions are not yet sufficient 
to completely determine all the neutrino relevant parameters. We will also briefly investigate 
the CP conserving three- zero textures, via both analytical and numerical analysis, to judge 
whether they are compatible with current neutrino data. 




(m ee ) 



i=l 



2 In the absence of CP violation, the quantities defined in Eq. ((HJ are modified to Xi = m, with i = 1, 2, 3, 4, 
and resultantly, the quantities in Eq. ([12]) are reduced to A;i = rrii/mi with i = 2, 3, 4. 
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III. ANALYTICAL ANALYSIS 



It seems impossible or unrealistic to obtain clean analytical expressions in the (3+1) 
model, as some of the elements of matrix U is rather lengthy. However, when we take into 
account the mass hierarchy m 4 3> mi : 2,3 and the smallness of the three ASM angles and # 13 , 
all the lengthy expressions can be sharply shortened and it becomes viable to analytically 
examine the 4x4 textures of M v with three zeros. The 120 textures has been classified 
and labelled in the Appendix. In this section we first analyze 99 of the 100 type I textures 
with zeros in the sterile sector and then the 20 of type II with all the three zeros in the 
active sector. For the texture compatible with experiments, we give explicit expressions 
and predicted ranges of the ASM angles and investigate the implications for the neutrino 
spectrum, CP violating phases, sum of neutrino masses and the effective neutrino mass in 
single and double beta decays. 



A. Textures of type I 

All the 100 type I textures with three zeros in the (3+1) model have sterile relevant 
vanishing elements, i.e., the vanishing elements in the fourth column and the fourth row. If 
the fourth diagonal element m ss of M v vanishes (class S), 

m 1 U 2 sl + m 2 U 2 s2 e ia + m 3 U 2 s3 e 1 ^ + m A U 2 s ^ = . (22) 

As 77ii ;2 ,3 7Ti4, to satisfy the equation, it is required that U si = cos#i 4 cos# 2 4 cos 6*34 <C 1, 
i.e., at least one of the three ASM angles be close to tt/2, which is not compatible with 



the small ASMs implied by SBL experiments, as mentioned in Ref. [25[. Furthermore, if 
#14 ~ 7r/2, it yields that 

mp ~ (m ee ) ~ m 4 ~ eV . (23) 
The current most stringent limit on (m ee ) is from the combined result of KamLAND-Zen 



and EXO-200, (m ee ) < (0.12 - 0.25) eV at 90% C.L. [34j, strongly disfavoring the value 
given in Eq. (|23|) . 

If one of the nondiagonal entries in the sterile sector is set to zero, i.e., m« = with 
i = e, fi, t (class J 7 ), 

rmUaUs! + m 2 U i2 U s2 e ia + m 3 U i3 U s3 e ip + m^JJ^ = . (24) 

For the same reason of 7711,2,3 <C m 4 , the equation requires Ui/JJ s ± < 1. As U e i ~ su, 
~ C14S24, U T 4 = C14C24S34, and U s $ = c u c 24: c 3 4, we arrive at the conclusion that, to 
satisfy Eq. (1241) , at least one of the two following conditions must be satisfied: (1) To have 
a vanishing entry m is (i = e, /1), the corresponding mixing angle between the active flavor 
i and the sterile one has to be nearly vanishing; or/and (2) At least one of the three ASM 
angles is close to it/ 2. Unfortunately, both of these two extreme conditions are excluded by 
SBL experiments. 

In the discussion above we did not consider the textures with solely m TS in the sterile 
sector, as, different from #14,24, a vanishing #34 is allowed due to the insensitivity of SBL 
experiments to v T . If 6*14,24 are mildly small, say ~ 10°, with the leading order (LO) approx- 
imation m 4 3> rni j2 ,3 and #1314 24 <^ 1, the third ASM angle is estimated from Eq. ff24|) to 
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be 



a/ I Am. 
\/ Am 



i 2 I 

734 7 =4^ ^4~0.4°. (25) 



In this case, the textures of class J 7 with m TS = cannot be completely excluded. However, 
such solutions are highly disfavored, as the parameter space in this specific region is rather 
small, which can be verified in the numerical calculation below. 

If the textures have both m ss = and m^ = with i — e, fi, r (classes A and C), as for 
class £, C14C24C34 <C 1 are required and all these texture are excluded by SBL experiments. 
In the texture of class B there is actually no mixing between the active and sterile sector 
and this texture is completely of no interest for us. For the textures of class T>, analogical 
to the class J 7 , two of the three ASM angles are required to be rather small or rather large, 
which are not allowed by SBL experiments. In short, all the (3+1) textures with zeros in 
the sterile sector are excluded by SBL experiments, except for the textures of class T with 
m TS = which, however, are highly disfavored. 



B. Textures of type II 

In three-neutrino scheme the M v textures with three zeros are definitely not compatible 



with the current neutrino oscillation data [18|. When the lepton sector of SM is extended 



to include one light sterile neutrino, one zero in the active sector of M v reads 

rrnUaUj! + m 2 U l2 U j2 e ia + m 3 U i3 U, 3 e i/3 + m 4 U l4 U j4 e^ = , (26) 

with i, j = e, /i, r. Compared to the three-neutrino framework, not only the fourth term is 
new, but the elements of U are tuned by six more neutrino parameters (three ASM angles 
and three more CP violating phases). All these facts provide the possibility of reconcile the 
experimental bound on neutrino parameters and the constraints from vanishing mass matrix 
elements. 

Numerically, 7711,2,3 <C 7774 and Uf 4 ~ 0.01 are generally much smaller than Ui m Uj n (i, j = 
e, fi, t, m, n = 1, 2, 3), thus it is possible for the contributions from active neutrinos to 
cancel out the fourth term in Eq. (1261) and obtain vanishing elements in the active sector of 
M v . If the sterile neutrino does mix with the active ones with strength of order 0.1 as the 
SBL experiments indicate, we can immediately obtain from Eq. (|2"6"|) that the active neutrino 
masses are of order m 4 L r l 2 4 ~ a/ Am^Uf^ ~ 0.01 eV and the lightest neutrino mass is of the 
same order 

m ~0.01eV. (27) 

Actually, this is a direct consequence of the neutrino mass matrix M v in presence of the 
sterile neutrino. To be more transparent, let us consider the simplified 2x2 version of M v 
with only the first generation of active neutrino, 

M v = ( ° mes V (28) 
V m es m ss J 

Due to the large active-sterile mass _splitting, we have set m ee = and are left with a mini- 



seesaw mechanism at the eV scale 35]. The ASM is expected to be ~ m es /m ss ~ 0.1; thus, 
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with an eV scale sterile neutrino m ss ~ 1 eV, we naturally obtain the active neutrinos are 
of mass ~ m 2 es /m ss ~ 0.01 eV. This meaningful prediction are not sensitive to the details of 
the active sector of M M (this part can always be ignored when we apply the seesaw formula) 
and it is a generic feature for 4x4 M v textures with zeros in the active sector, regardless of 
the number of zeros {25|. One phenomenological consequence of the prediction on neutrino 
spectrum is that models yielding a massless state or large mass hierarchy in the active 
neutrino sector are not compatible with the (3+1) scenarios with texture zeros. 

Applying the approximations mi i2 ,3 <C m 4 and #13,14,24,34 <C 1, the six independent 
vanishing elements in the active sector of M v are simplified to, at LO, 
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(29a) 
(29b) 
(29c) 
(29d) 
(29e) 
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We emphasize again that these explicit equations strongly indicate the feasibility of the 
(3+1) model with zeros in the active sector: Assuming mo ~ 0.01 eV, Eq. (1291) implies that 
the ASM angles 

1/2 

-(0.1-0.2), (i = l, 2, 3) (30) 



SiA 



V Am 21 


1/2 






or 





coinciding with the hints of SBL experiments. In addition, a few comments on the approxi- 
mate expressions follow: 

• U e 4 ~ S14, thus after approximation the vanishing condition m ee = is reduced to the 
LO expression for s 2 4 , i.e., (I29ap . Similarly, by means of ~ s 2 4 and U r4: ~ s 34 , 
the two vanishing diagonal elements = m TT = are simplified to the equations 
for the angles #24,34 [( I29bl) and ( j29cf) ] . while the three nondiagonal vanishing elements 
m efl = m er = m^r = lead to the three equations ( I29dl) . (I29el) . and ( I29fl) of the cross 
terms. 



In Eqs. ( I29dl) and (I29el) the terms of form 7723513 are reserved, because for the neutrino 
spectrum of NH their magnitudes are of the same order of the m 12 terms, whereas in 
the other four equations of (1291) . the LO terms in #13 are of form s 2 3 and are neglected. 

Comparing the Eqs. f!29bl) and f|29cp as well as f!29dl) and f|29ep . these equations exhibit 



the approximate fi — r symmetry [23[ . For example, the textures Mv with m 
can be approximately obtained by acting the fi — r permutation matrix 

/o o\ 

10 



P. 



10 




(31) 



12 



on the matrices Mu^ with m 
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MP = PrrMV>P*. 



(32) 



This approximate symmetry will be more transparent if we write explicitly down the 
analytical expressions of the ASM angles, as shown below. 

All the discussions of phenomenological implications of the 20 type II texture are based 
on the equations in (|29|) . With the six equations, we can easily obtain clear expressions of 
the quantities A21, 31,41 and then immediately the neutrino mass ratios and the Majorana 
phases. However, in the (3+1) model it is more attractive to investigate directly the ASM 
angles and relevant phenomenologies. We denote the six vanishing conditions of Eq. ( 12 9 p 
simply as 

(33a) 
(33b) 
(33c) 
(33d) 
(33e) 
(33f) 



where 
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Sl4 


= a , 




= 6, 


■S34 


= c, 


S14S24 


= d 2 


S14S34 


= e 2 


524^34 


= f 



1 



m 4 
1 



'm 4 



24, 24,o 22 1 

m 1 c 12 + m 2 s 12 + 2mim 2 s 12 c 12 cosaj 



1/4 



244, 244, 24,0 224 
™1 S 12 C 23 + m 2 C 12 C 23 + m 3 S 23 + 2m 1 m 2 S 12 C 12 C 23 COS Q 

+2mim 3 s^ 2 s 23 c 2 3 cos (3 + 2m 2 m3C^ 2 s 23 c 23 cos(a — /?)] 



1/4 



,2„4 4 



2„4 4 



,2„4 



" m l s 12 S 23 ~l~ m ') c W s ")^ + m S c 



2 U 12°23 



'3^23 



2mim 2 s^ 2 c^ 2 s 23 cos a; 



1/4 



1/4 



+2m 1 m 3 s 12 s 23 c 23 cos /3 + 2m 2 m 3 c 12 s 23 c 23 cos(a 

/ 2 , 2\ 2 2 2 , 222 <-> 222 

[m x + m 2 )s 12 c 12 c 23 + m 3 s 23 s 13 - 2m 1 m 2 s 12 c 12 c 23 cos a 
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(34a) 



(34b) 



(34c) 



(34d) 



(34e) 



(34f) 



The expressions in Eq. ( 134"|) and then the explicit formulas for the ASM angles ( 1331) can be 
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viewed as expressions of the lightest neutrino mass mo, with, for the NH spectrum 



nil — mo 



m 3 = s jml + Am^ , 
WI4 = \j TOq + Am 



2 

41 • 



or for IH 



nil = \/ m o + A m !i , 



m 2 = ymg + Am^ + Am^ , 

m 3 = 

ni4 = \J wig + Am| x + Amli . 



(36) 



In the case of IH, as m 3 = m ~ 0.01 eV and consequently the terms m3Si3 are much smaller 
than the contribution from m^, the d and e in Eq. (134|) can be further simplified to 

d = ——= [(ml + ml)sl 2 cl 2 cl 3 — 2niim 2 s\ 2 c 2 l 2 c\ z cosa ] ± ^ > (37a) 
\i 7724 



[(m\ + ml)sl 2 c 2 i 2 sl 3 - 2mim 2 sl 2 c 2 i 2 s 2 23 cos a] 1/4 . (37b) 



With the four analytical expressions f l34b|) - fl34e]) for the ASM angles ■524; S34, S14S24, and 
S14S34, the approximate \i — r symmetry is clearer than in Eq. ff29l . 



#12 - > #12 , #13 813 , #23 - > f — #12 
a — > a , j3 — > (3 , <5i3 — ^ — tt . 



(38) 



similar to the case for the schemes with three flavors of neutrinos [23]. It should be reminded 
that the simple transitions fl38|) are obtained based on the assumptions of small reactor 
mixing and ASMs #13,14,24,34 <C 1; the exact formula is a bit more complex, as presented in 



Ref . [25J . For the 20 type II textures we now begin to analytically analyze, the approximate 
fi — r symmetry relates the textures Bi — B 3 , B 2 — B 4 , B 5 — B 6 , Ci — C2, Di — D 2 , D 3 — D 5 , 
D 4 — D 6 , and E 2 — E 3 . In our analytical analysis, we consider all the 20 textures to check 
whether the approximate symmetry holds, as we expect, for the eight pairs of textures. 

• Texture A 

The first three diagonal elements of M v vanishes and this texture involve the first 
three expressions of the ASM angles in Eq. (|33|) . Assuming niQ = 0.01 eV, with the 
3(7 ranges for the neutrino mass squared differences Am^ 31 41 and th mixing angles 
#12,23,13 given in Tables HI and HT1 as well as a, (3 G [0, 2n], one can easily calculate the 
allowed ranges of the three ASM angles for both the two neutrino spectrums of NH 
and IH, 

1.6° < #14 < 6.6° (NH) or 5.0° < 6 M < 14.4° (IH) , 

4.5° < # 24 < 12.6° (NH) or < # 24 < 12.0° (IH) , (39) 

3.9° < #34 < 12.3° (NH) or < # 34 < 12.3° (IH) . 

The predicted ranges substantially overlap with the ranges determined by SBL exper- 
iments, which reveals that both the NH and IH spectrums are allowed in the texture 
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A. The similar predictions on the angles #24,34 indicates that A is approximately self- 
symmetric under the \i — r permutation P Atr . The three textures C3, Ei, and F are 
also self-symmetric under the approximate fi — r permutation, as shown in Tables II III 
and IIVI below. In Eq. (1391) . the lower limits of the angles #24,34 for IH vanish, because 
in Eqs. (134bj) and (I34cl) it is possible for the contributions from the three light states 
m i,2,3 to completely cancel out each other. In addition, Eqs. f )34ap — f)34cp imply that 
it is easier to produce vanishing elements (m eejWJ TT = 0) if the Majorana CP violat- 
ing phase a ~ ir than a ~ 0, especially for the spectrum of IH. The reason is that 
with smaller values of the ASM angles, it is more likely for the light states in 
Eq. fl29l) to cancel out the contribution from the sterile state 7714. Although the Majo- 
rana phases cannot be observed in neutrino oscillation experiments, they are relevant 



to the effective neutrino mass in Qvf3f3 [36|. We now skip over the discussions of the 



implications on the effective neutrino mass in single and double beta decays; at the 
end of this subsection we will collectively analyze these phenomenological implications 
for all the type II textures. 

Class B 

Let us first take the texture Bi as an explicit example, in which m eT = in addition to 
m ee = m MA1 = 0. The LO expressions for S14 and S24 are the same as for texture A, and 
s 13 can be solved from Eq. (I34ep . s 34 = e 2 /si 4 = e 2 /a. Assuming again m = 0.01 eV 
as done for texture A, we can immediately obtain the predicted ranges of the ASM 
angles: 

1.6° < #14 < 6.6° (NH) or 5.0° < 9 14 < 14.4° (IH) , 

4.5° < #24 < 12.6° (NH) or < # 24 < 12.0° (IH) , (40) 

< #34 < 26.6° (NH) or < # 34 < 23.7° (IH) , 

which are also compatible with the SBL constraints. Due to the minus sign in Eq. ( 137bj) 



for the IH spectrum, when the Majoran phase a ~ it, the angle #34 tends to take larger 
values than with the assumptions of a ~ (although the other ASM angles #14, 24 are 
likely to be smaller), thus the IH spectrum is less favored than NH. 3 Due to the same 
reason of preference for small ASMs, the Majorana phase a ~ ir is much more favored 
than a ~ 0. 

The analysis above apply analogically to the textures B 2 — B 6 . The compact analyt- 
ical expressions and the resultant numerical ranges of the three ASM angles in these 
textures are listed, respectively, in Tables HTT1 and HVl In both tables, the textures with 
the same superscripts are related by the approximate \i — r symmetry, under which 
the analytical expressions simply transform as 

b< — >c, d i — >e. (41) 

The texture B 3 is the partner of B x under the \i — r symmetry. Comparing their 
analytical formula and numerical predictions in Tables IIHI and IIV[ it is clear that the 
fi — t symmetry approximately holds, as we expect. These two textures share many 
similar features: B 3 also prefers much for NH and a ~ ir. 



As explained for the texture A, with larger values of 634, it is harder to produce a vanishing m er . Although 
the SBL constraint on 634 is less severe than ^14, 24, it is completely irrelevant to the discussion here. 
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TABLE III: LO expressions of the three ASM angles for the 20 textures of type II. The angle 
#34 in the texture Ci, #24 in C 2 , and #14 in C3 are unconstrained. The textures with the same 
superscripts are related by the approximate [i — r symmetry. See text for details. 
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Similarly, it is easy to check that the approximate [i — r symmetry holds true for the 
two other pairs, i.e., B 2 — B4 and B5 — B6, and we can predict that both NH and IH 
are compatible with the four structures. With regard to the CP violating phases, the 
textures B 2j 4 prefer the values of a around ir. When it comes to Bs^, the situation 
becomes a bit more complex. In the two textures, both the two vanishing elements 
m p,p,,TT — involve the two Majorana phases a and /3 at LO, as shown in Eq. (134bj) 
and (I34cl) . Among the four approximations of the two phases 

a~0&/3~0, a~7r&/3~7r, a~0&/3~7r, a~7r&/3~0, (42) 

for the NH spectrum, small ASMs are most acceptable when the third approxima- 
tions is satisfied, in which case the contributions from the two states mi i2 interfere 
destructively with the third state m 3 . For IH, we do not have such clear predictions. 

• Class C 

We first consider the texture Ci, in which m efJL = in addition to m ee = m MAt = 0. 
The three corresponding Eqs. fl34al) . f!34bl) . and f!34dl) are irrelevant to the ASM angle 
6* 3 4; actually even the three exact vanishing conditions m ee>flfl)TT = without any 
approximation does not involve # 34 and thus #34 are completely unconstrained. The 
predicted ranges of 6 , 14j24 are exactly the same as in the texture A, assuming m = 
0.01 eV, 

1.6° < 0i4 < 6.6° (NH) or 
4.5° < 9 24 < 12.6° (NH) or 

whereas Eq. f ]34dj) implies 

< v/rai < 7.6° (NH) or 

Obviously, these three vanishing conditions are consistent with each other. When 
the neutrino spectrum is inverted, the vanishing nondiagonal element m efJ- is largely 



5.0° < 9 14 < 14.4° (IH) , 

< 6 24 < 12.0° (IH) , 1 ' 

< V9^b\i < 12.8° (IH) . (44) 
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TABLE IV: Numerical ranges of the LO analytical expressions given in Table IIIII for the ASM 
angles in the 20 textures of type II. "— " denotes that the corresponding angle is unconstrained. 
The textures with the same superscripts are related by the approximate fi — r symmetry. See text 
for details. 
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[0, 26.6°] 
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incompatible with the two diagonal conditions: When a ~ ir Eq. (I37ajl produces larger 
values of 9^ 2 4 than Eqs. (I34aj) and ( I34bl) . whereas when a ~ 0, it leads to smaller 



values of #14, 2 4. Thus IH is highly disfavored by the texture Ci. For the NH spectrum, 
however, the three vanishing conditions are compatible and the phase a ~ ir is favored 
due to the same reason as stated for the textures A and Bi. 

At LO, the texture C2 does not involve the angle 624 and C3 is irrelevant to 614. 4 The 
expressions and numerical ranges of the other ASM angles for C2 and C3 are given, 
respectively, in Tables DTD and [TV] C2 is related to Ci by the /i — r symmetry, and 
it favors the NH spectrum and the values of a around ir. For texture C3, as a result 
of the different signs of the cross terms in Eq. (I34fj) from that in (I34bj) and ( I34c[) . no 



4 



At next-to-leading (NLO) order, the two textures C2, 3 involve all the three ASM angles, but the formulas 
are rather complex. Thus we do not show the formulas of S24 of C2 and S14 of Ci in Table IIIII The 
preferred ranges of the two angles can be easily recognized in the numerical calculation below. 
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matter (3 is ~ or ~ n, it is challenging for ( I34fj) to keep consistent with ( 134b|) and 
f)34cp and then this texture is highly disfavored. 

Class D 

In the first texture Di of this class, m ee = and m efl = = 0. We can easily obtain 

the analytical formulas of the three ASM angles from Eqs. f)33ap . (133d[) . and ( I33fj) . 



S14 = a , (45a) 
s 2 4 = d 2 /a , (45b) 
s 34 = af/d 2 . (45c) 

Assuming mo = 0.01 eV, the numerical ranges of th ASM angles are, 

1.6° < 0i4 < 6.6° (NH) or 5.0° < 9 U < 14.4° (IH) , 

0.7° < fl 24 < 26.6° (NH) or 0.7° < # 24 < 23.0° (IH) , (46) 

2.4° < fl 34 (NH) or 0.2° < # 34 (IH) . 

Originally, the lower bounds on # 2 4 for both NH and IH spectrums are vanishing. 
However, as the constraint on #34 is obtained by means of S34 oc l/s 2 4, to ensure that 
S34 < 1, a cutoff is imposed on the lower end of s 2 4, which is estimated to be 0.7° as 
shown in Eq. (146|) . Due to the opposite sign of the cross term mim 2 in Eqs. (I34ap 
and (I37al) . the IH spectrum is highly disfavored by the texture Di; for NH the phase 
a ~ 7r is much more favored than a ~ 0. 

For the textures D 2 — D 6 , see Tables DTD and [IV] for the expressions and expected 
numerical ranges of the ASM angles, which reveal the relation of Di — D 2 , D3 — D5, 
and D4 — D6 under the approximate fi — r symmetry as we expect. With the method 
of rough estimation frequently applied above, we can foresee that only the texture D 2 
disfavor IH. In addition, the textures D 2> 4> 6 require that a ~ n, while D 3i 5 are more 
acceptable if a ~ and (3 ~ w. 

Class E 

In the texture Ei, m ee = and m efl = m eT = 0. After applying Eq. (I33aj) . s 2 4 and S34 
can be written as functions of Su and then 

S14 = a , (47a) 
s 24 = d 2 /a, (47b) 
S34 = e 2 /a . (47c) 

With the lightest neutrino mass mo = 0.01 eV, estimated ranges of the ASM angles 
are, for both NH and IH spectrums, 

1.6° < #14 < 6.6° (NH) or 5.0° < 6 U < 14.4° (IH) , 

< #24 < 26.6° (NH) or < # 24 < 23.0° (IH) , (48) 

< #34 < 26.6° (NH) or < # 34 < 23.7° (IH) . 

Based on the argument analogical to that for the texture Di, IH is high suppressed 
compared to NH, in which the preferred value of a is around ir. 
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The expressions and numerical ranges of the ASM angles in the textures E2, 3 are pre- 
sented, respectively, in Tables HTT1 and HVl Obviously, the approximate fi — r symmetry 
holds for the two textures, and both the two neutrino mass hierarchies are compatible 
with the two structures. For the NH spectrum, to obtain small ASM angles, a ~ 
and /3 ~ 7r is required in both the two textures, while for IH, a tends to be around tt. 

Texture F 

In this texture the three nondiagonal elements in the active sector of M v vanish, i.e., 
m efJi = m er = m^ T = 0. Based on Eqs. ( 133dj) — (133fj) . the three ASM angles can be 



expressed as 

s u = de/f, 
S24 = df/e, 

s M = ef/d, (49) 

Assuming as above m = 0.01 eV, we can estimate the ASM angles, 

0.3° < #14 < 4.4° (NH) or 0.1° < 6> 14 < 70.8° (IH) , 

0.8° < # 24 (NH) or 0.9° < 6 2i (IH) , (50) 

0.7° < #34 (NH) or 1.0° < 6 3i (IH) . 



As seen in Eqs. fl34dj) -f l34f|l . all the three CP violating phases a, 0, and #13 play an 



important role in pinning down the neutrino spectrum. It is not so transparent as the 
textures above to figure out the preferred neutrino spectrum and the favored regions 
of the phases. More phenomenological discussions of this textures are postponed to 
the next section, wherein the numerical analysis can help us to further investigate this 
texture. 

After the detailed investigation of the neutrino spectrum, ASM angles and CP violating 
phases of the 20 type II textures, we now estimate both the effective neutrino masses in 
both tritium beta decay and Ovfift. With mo ~ 0.01 eV, the contribution from the sterile 
state dominates the beta decay electron neutrino mass 



m p = yj miclzcizcii + m|s? 2 c? 3 c?4 + m|sf 4 ~ 0.1 eV . (51) 

In Eq. ( l5"Tj) we neglect the contribution from m 3 , as its is highly suppressed by the smallness 
of reactor mixing sf 3 . Assuming explicitly mo = 0.01 eV as well as the 3a ranges of the active 
neutrino parameters and Am 41 and # 14 listed in Tables [J and [Til one can more accurately 
predict the mass mp in general (3+1) neutrino models: 

91 meV < m p < 420 meV (NH) or 102 meV < mp < 423 meV (NH) . (52) 

With the help of Eqs (1331) and ( 1341 . we calculate the effective mass mp in all the 20 textures, 
as collected in Table |V] For both NH and IH spectrums, mp in all the textures are of order 
0.1 eV, consistent with the rough estimation. Furthermore, the predictions presented in 
the table overlap partially with Eq. (152]) . which implies again that these type II textures 
are compatible with experiments. The lower limits ~ 13 meV and ~ 48 meV correspond 
to a vanishing 6 , 14 and denote the minimum of the contribution from the active neutrinos 
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TABLE V: LO estimations of the effective electron neutrino mass mp in beta decay (left table) 
and the effective Majorana electron neutrino mass (m ee ) in Ov/3/3 (right table), with all in units of 
meV. See text for details. 



Textures 


NH 


IH 


Textures 


IN xl 1x1 


A, Bi 2,3,4, Ci 2, Dl,2j Ei 


36 - 135 112 - 294 


A, Bi ,2,3,4, C12, Di,2j Ei 





B5, D3, E2 


13 - 152 


>48 


B5, D3, E2 


- 15 > 47 


B6, D5, E3 


13 - 166 


>48 


B6, D5, E3 


- 16 > 47 


D 4 


13- 166 


>48 


D 4 


- 22 > 47 


D 6 


13- 159 


>48 


D 6 


- 21 > 47 


F 


13 - 90 


>48 
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- 13 > 47 
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with the NH and IH spectrums, respectively. If the neutrino spectrum is inverted, rxip is 
mildly augmented as a result of heavier 77*1,2- The ASM angle #14 in the texture C3 are not 
restricted at LO, thus the neutrino mass mp for C3 are left blank in Table IVl The designed 



discovery sensitivity of the KATRIN experiment is about 0.35 eV [37[ and it seems possible 
for KATRIN to distinguish the neutrino mass hierarchies in these textures, if a positive 
signal is observed in this experiment. 

With regard to the effective Majorana neutrino mass (m ee ) in 0^/3/3, if mo = 0.01 eV, 
neglecting the highly suppressed m 3 s\ 3 term, 

(m ee ) = m ee = \rriic 2 12 cl 3 c 2 u + m 2 s\ 2 c\ 3 c\ i e ta + m 4 Si 4 e i( " 2<5l4+7) | . (53) 

Rough estimation of the order of magnitude tells us that (m ee ) is of order 0(10 me V), which 
is just of the same order as the sensitivity of the next generation Ov/3/3 experiments [36] . The 
leading order predicted ranges of (m ee ) for the 20 textures are listed in Table |V] The ten 
textures A, Bi 2,3,4, Ci,2, Di 2 , Ei assume m ee = 0, thus definitely they cannot produce 
any signal in the 0z//?/3 experiments. With help of the Majorana phases a and 7, the nine 
textures B 5 6 , D 3i 4 j5i6 , E 2 ,3, F with NH spectrum can lead to vanishing Majorana neutrino 
mass, and the maximum value of (m ee ) in these textures at the leading order is about (10- 
20) meV. For IH, m 12 become heavier and it is impossible for the states m 12i4 to completely 
cancel out each other, thus we arrive at a nonvanishing lower bound of ~ 47 meV in the 
nine textures. At LO the texture C3 cannot predict (m ee ), however, if it is compatible 
with the SBL experiments, the prediction of C3 on (m ee ) is about the same as the nine 
textures such as B5. The most meaningful phenomenological features in Table M is that 
with m = 0.01 eV the NH spectrums of the ten textures B 5 6 , C3, D 3i 4 i5i 6, E 2 ,3, F are 
clearly distinguishable from the IH spectrums, and it is promising to test the neutrino mass 
hierarchies in these textures by means of the next generation Qvflfl experiments [36} . The 
phenomenological implications of these textures on the single and double beta decays will 
be further investigated in the following section on numerical analysis. 



IV. NUMERICAL ANALYSIS 

We now turn to numerical analysis of the 120 textures of M v with three zeros in the 
(3+1) neutrino model. In this section we first state in detail how the numerical calculation 
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is implemented and then analyze systematically the numerical results of the type I and 
II textures. A number of scatter plots are depicted to illustrate the phenomenological 
consequences of these textures: neutrino mass spectrum, ASM angles, CP violating angles, 
and the nonoscillation quantities, i.e., sum of neutrino mass E m , effective neutrino mass mp 
in beta decay and the Majorana neutrino mass (m ee ) in 0u(3(3. 

A. Calculation procedure 

The numerical calculation is implemented as follows: 

• The three known active neutrino mixing angles #12, 23, 13 are allowed to vary randomly 
in their 3a ranges presented in Tabled! the three ASM angles #14,24,34 and three Dirac- 
type CP violating phases $13, 14, 24 are viewed as "unknown" parameters and they vary 
randomly in the range [0, 7r/2] and [0, 2tt], respectively. 

• For each randomly generated group of the nine physical parameters given above, we 
calculate the corresponding A2i,3i,4i and the two resultant neutrino mass square dif- 
ference ratios R v \ and R v2 , based on Eqs. (EJ), (TT21) . (TT3"j) . and f[T7l) . At this stage, the 
three neutrino mass ratios m 2 /mi, 777.3/mi, m 4 /mi and three Majorana-type phases a, 
(3, 7 can be obtained according to, respectively, Eqs. ( ITS]) and (TIB]) . If R vl and R u2 fall 
simultaneously inside the ranges determined by the data on mass square differences in 
Tables H and M 

0.0260 < R ul < 0.0371 , 

2.81 x 10~ 5 < R u2 < 1.14 x 10" 4 , (54) 

then the corresponding input groups of neutrino parameters are possibly compatible 
with experimental constraints, otherwise, they are excluded by experiments at 3a C.L. 

• For the groups of neutrino parameters consistent with experiments obtained above, 
we randomly generate Am^ in its 3a range. With the mass ratio m 2 jra\ obtained 
above, we can easily get the corresponding absolute neutrino masses m 12 , and further 
produce m 3 4 when the mass ratios m 3 /mi and m 4 /mi are also involved. Given the 
four neutrino masses ?n 12 ,3,4, it is compulsory to check if the corresponding jAm^l 
and Am^j do fall into their corresponding 3a ranges and discard those that do not. 

• Finally, we obtain all the 16 neutrino physical parameters, among which six are con- 
strained by experimental data: 6*12,23,13, Am^, jAm^l, and Am| 1; six are randomly 
generated unconstrained: #14,24,34 and #13, 14, 24, and the last four, i.e., a, (3, 7, and 
mo, are obtained based on the correlation equations implied by the three vanishing 
elements of M v . 

As aforementioned, we own 10 "unknown" parameters and six constraints, it seems more 
smart and economical if we have four, instead of six, unconstrained parameters in our 
calculation. However, it is unclear how to choose four from the ten "unknown", and more 
important, obtaining the six "unknown" parameters from the correlation equations ( T151) and 
f fT6|) would become a troublesome task. Comparatively, our numerical calculation is more 
straightforward and feasible. Our calculation can be regarded as a method in which two of 
the six correlation equations are solved with the Monte-Carlo method [jjj]. 
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As seen in Table [J the 3a range of #23 for NH spectrum is slightly different from the case 
of IH, and the difference of 3a ranges of lAmf-J for the two mass hierarchies is somewhat 
even larger. Without bias on any of the two hierarchies, at the beginning of numerical 
calculation we set the covering ranges for both NH and IH as the 3a ranges of jAm^l and 623- 
Numerically, they are (2.21 — 2.74) x 10~ 5 eV 2 and 0.36 — 0.68, respectively. After obtaining 
the neutrino masses by means of the calculation procedure above, we can immediately 
identify the mass hierarchy for each group of neutrino parameters. We check whether the 
mass hierarchy and the values of jAm^l and 823 in a given group of neutrino parameters 
are allowed by the corresponding experimental data listed in Table HI and discard those that 
are not. 



B. Textures of type I 



As explained in the analytical analysis, the textures with zeros in the sterile sector require 
large or nearly vanishing ASM angles due to the large mass hierarchy between the active 
and sterile neutrinos <C m^. We numerically calculate all the 100 type I textures and 

find that all the significant features are consistent with the analytical expectations. 

As a concrete demonstration, the randomly generated points of the ASM angles of the 
two textures 
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(55) 



are plotted in Fig. [TJ In the texture a with m ss = 0, it is transparent that the ASM angle 
#14 is very large or, less likely, #24 is close to the maximum value 7r/2, consistent with the 
analytical expectations j25[. For the texture b, two nondiagonal elements in the sterile sector 
are zero, i.e., m es jTS = 0. We can easily read from the middle and right panels of Fig. [T]that 
#14 is mostly very small and the other two ASM angles are generally very large, which is just 
predictions of the analytical analysis. As stated above, the SBL experiments are insensitive 
to the tauon flavor and thus permit a vanishing #34. The matrix 
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(56) 



has only one zero m TS = in the sterile texture, and the numerical predictions of the three 
ASM angles are shown in Fig. [2j As we expect, #34 is mostly close to 0, while 624 is mostly 
close to 90°. By examining the neutrino parameters produced in the numerical calculation, 
it turns out that only rarely can we obtain the three ASM angles all compatible with SBL 
experiments, i.e., this texture is highly disfavored, although it has not been completely 
excluded. 



22 



abb 




014 [deg] M [deg] <9 24 [deg] 



FIG. 1: Distributions of the ASM angles for the textures a and b given in Eq. (|55H . The left panel 
is the #14 — 624 distribution in the texture a, in units of degree; the middle and right panels show, 
respectively, the 6*14 — #34 and 824 — #34 distributions for the texture b. NH and IH are denoted, 
respectively, as blue filled circles and red squares. 




014 [deg] 0,4 [deg] 

FIG. 2: Distributions of the ASM angles for the texture c given in Eq. (I56p . The left and right 
panels are, respectively, 614 — 624 and 614 — 634 distributions, in units of degree. NH and IH are 
denoted, respectively, as blue filled circles and red squares. 

C. Textures of type II 

The analytical analysis reveals that, without compulsorily requiring nearly vanishing or 
very large ASM angles, it is possible for the type II textures to be compatible with SBL 
experiments. Now we present the numerical results of the 20 textures, by sort of the various 
phenomenological implications: 

• The most significant feature which are common to all the 20 textures is that the 
lightest neutrino mass m for both NH and IH is of order (0.001 — 0.1) eV, as clearly 
shown in Fig. |3j This coincides with the analytical prediction that m are required to 
be around the order 0.01 eV. 

• In Fig. [31 we present the numerical predictions of the three ASM angles of texture A 
for both NH and IH spectrums. It is transparent that the ranges of the ASM angles 
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A : NH A : IH 




m [eV] m [eV] 

FIG. 3: The red filled circles, green squares, and blue diamonds denote, respectively, #14, 624, and 
634 in the texture A as functions of the lightest neutrino mass mo, with NH (left panel) or IH 
(right panel) spectrum. 



with mo = 0.01 eV are almost exactly the same as the analytical predictions. On the 
whole, if m Q is larger, the ASM angles tend to be larger, in contrast, the angles are 
smaller for smaller values of m . The figures of the ASM angles for the textures B 1; 
Ci, Dx, Ex, and F are shown in Fig. HJ It is possible that NLO and higher order 
corrections beyond Eq. (129]) can widen somewhat the LO allowed ranges of the ASM 
angles. For instance, in the upper left panel of Fig. [4j when mo = 0.01 eV, the angle 
834 of Bx with NH spectrum can take values larger than the LO upper bound 26.6° 
given in Table IIV1 However, the simulations are consistent with the analytical analysis 
on the whole, which not only verifies the feasibility of the LO analytical analysis, but 
demonstrate again that the type II textures are more or less compatible with the active 
and sterile neutrino data. 

The three vanishing elements of the texture Cx do not involve 8 34 , thus this angle is 
completely constrained. In our simulation it distributes evenly in the range [0, vr/2], as 
we expect, shown in the upper left panel of Fig. The other three panels of Fig. are 
on the angle 824 of texture C2 (NH) and 614 of C3 (NH and IH) which are restricted 
at NLO. The numerical calculation reveals that 824 show a slight preference for larger 
values around 80°, whereas 814 are mostly of the values 70° — 90°. Consequently, the 
texture C3 is excluded by SBL experiments. In all the discussions below, we will not 
further consider the texture C 3 . 

• With regard to the neutrino spectrums of the 19 textures, the numerical calculations 
are also consistent with the analytical predictions: the NH spectrum is much more 
favored by the seven textures Bx,3, Cx,2, Di 2, Ex than IH, whereas both NH and IH 
are allowed in the other twelve textures A, 62,4,5,6; D 3 4 5 6 , E 2 ,3, and F. We delib- 
erately point out the mass spectrum of F cannot be clearly determined by analytical 
analysis. 

• Analytical analysis reveals that some of the CP violating phases are severely con- 
strained, especially the Majorana phases a and 0. The preferred ranges of the Dirac 
and Majorana phases obtained in the numerical calculation are collected in Table IVII 
This table tells us that not only the feasibility of the analytical analysis as well as the 
predictive power of the type II textures on the phases are verified, but more predictions 
can be yielded in the numerical calculation on the phases, i.e., the predictions of #13,24 
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FIG. 4: The red filled circles, green squares, and blue diamonds denote, respectively, #14, #24 > and 
#34 as functions of the lightest neutrino mass too, in the textures Bi (NH, upper left panel), Ci 
(NH, upper right panel), Di (NH, middle left panel), Ei (NH, middle right panel), and F (lower 
left panel for NH and low right panel for IH). 



and 7 presented in the table. In addition to the Dirac phases #13, the two Majorana 
phases a and (3 in the texture F cannot be clearly predicted by analytical analysis, 
as we have mentioned above. As representative examples, we present in Fig. the 
numerical results of the phases #13 in the texture F (NH), 624 in B2 (NH), a in A 
(IH) and D 5 (NH), (3 in E 2 (NH), and 7 in D 5 (NH). All the numerical results are 
explicitly consistent with the analytical predictions. 

• As the analytical analysis, the numerical calculation can also predict the physical 
nonoscillation observables: sum of neutrino masses S m , the effective neutrino mass mp 
in single beta decay and the neutrino mass (m ee ) in Oz//3/3. As an explicit illustration, 
the three quantities for the texture D3 are presented in Fig. [7J The three light active 
neutrinos are generally of order 0.01 eV and the sterile one of order 1 eV, thus the 
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FIG. 5: Distributions of the ASM angle #34 of the texture Ci (NH, upper left panel), #24 of C2 
(NH, upper right panel) and #14 of the texture C2 (lower left panel for NH and lower right panel 
for IH) as functions of the lightest neutrino mass niQ. The blue and red points denote, respectively, 
the NH and IH spectrums. 



TABLE VI: Preferred values of the CP violating phases in the 19 type II textures obtained in the 
numerical calculation. 



Phase Preferred values Textures 


613 tt/2 or 3vr/2 


F 




<5 24 


B2,4, ^4,5,6! ^2, 


3 and F, all with NH spectrums 


or 7T 


A (NH); F (IH) 




vr/2 or 3tt/2 


Ei 
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D 3 , 5 ; B-5,6 (NH), 


E 2 , 3 (NH) 


7T 


A, Bi ,2,3,4, C12 


i Di,2,4,6, Ei, F; E2,3 (IH) 
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D 3 , 5 ; B 5 ,6 (NH), 


E 2 , 3 (NH) 


7 7T 


D 5 





masses of the four states sum up to ~1 eV, as shown in the upper panel of Fig. [7J which 
is compatible with the bound on neutrino mass from cosmological observations |||. 

Numerical prediction of the effective neutrino mass nip in tritium beta decay for the 
texture D3 is presented in the lower left panel of Fig. [7J Considering the SBL con- 
straints, the lightest neutrino mass mo cannot too larger than 0.01 eV, or the ASMs 
would become very large. Thus, assuming m — 0.01 eV, the effective neutrino mass 
mp is around (^(lO^eV), consistent with the analytical prediction in Sec. IIHI A more 
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FIG. 6: Distributions of the CP violating phases <5i3 in the texture F (NH, upper left panel), 624 
in B2 (NH, upper right panel), a in A (IH, middle left panel) and D5 (NH, middle right panel), 
f3 in E2 (NH, lower left panel) and 7 in D5 (NH, lower right panel), as functions of the lightest 
neutrino mass ttiq. The blue and red points denote, respectively, the NH and IH spectrums. 



phenomeno logically meaningful output is that the upper bound on mp for the NH 



spectrum is below the KATRIN discover sensitivity of 0.35 eV [37j, as indicated by 



the LO analytical expressions. If the electron neutrino mass mp is above 0.35 eV, 
KATRIN can in principle distinguish the two neutrino spectrums. However, the pa- 
rameter space for the NH spectrum with mo — 0.01 eV and 0.35 eV < < 0.423 eV 
is rather small and it is challenging for KARTRIN to produce positive signal in this 
restricted region and finally pin down the neutrino spectrum. If the KATRIN experi- 
ment cannot detect electron neutrino mass but otherwise set a upper limit of 0.2 eV, 
then we are not able to differentiate the two neutrino hierarchies. 

For Ovfifi, as shown in the lower right panel of Fig. [7J with m ~ 0.01 eV, the effective 
Majorana neutrino mass (m ee ) for IH spectrum is large than ~ 40 meV, consistent 
with the analytical predictions. The (m ee ) for NH spectrum is clearly much lower, 
predictively around 10 meV. The manifestation of nonvanishing lower limit for NH, 
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FIG. 7: Distributions of the physical nonoscillation quantities for the texture D3 as functions of the 
lightest neutrino mass mo- sum of neutrino masses S m (upper panel), effective neutrino mass mp 
(lower left panel) in beta decay and the effective Majorana neutrino mass {m ee ) (lower right panel) 
in Ov/3/3. The blue filled circle and red square denote, respectively, the NH and IH spectrums. 



which seems not consistent with the LO prediction, is a result from the corrections of 
the NLO and higher order terms. As mentioned in Sec. IIII[ its is promising for the 
next generation 0v/3{3 experiments to pin down the neutrino mass spectrum and test 
the type II textures jltj ]. 

The LO analytical analysis indicates that there exist approximate \l — t symmetry between 
some of the type II textures. In the numerical calculation without any approximation, the 
approximate \l — t symmetry are reserved and the textures related by this flavor symmetry 
share the same phenomenological features, e.g., the preference of the textures Bi 3 for the 
NH spectrum over IH. One antiexample is that the texture D5 prefer the CP violating phases 
#24 ~ and 7 ~ Ti but its partner D 3 does not, as presented in Table fVTl We speculate this 
explicit and significant flavor "asymmetry" is due to the NLO and higher order corrections 
beyond Eq. fl22). 



V. TEXTURES WITHOUT CP VIOLATION 

When all the CP violating phases vanish, i.e., #13,14,24 = a = /3 = 7 = 0, without tuning 
of the phases it would be more challenging for the vanishing conditions ( TTUj) to be satisfied. 
For instance, the CP conserving partners of the six vanishing conditions in the active sector 
of M v are, respectively, 

m ee : m 4 s^ 4 + ( mic^ + rn 2 sl 2 ) = , (57a) 
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: m 4 S24 + [rriis{ 2 (% 3 + m 2 c\ 2 c 2 23 + m 3 S2 3 J = , (57b) 

m TT : m 4 S34 + \m\s\ 2 s\ 3 + m 2 Ci 2 S2 3 + m 3 cl 3 j = , (57c) 

m e/i : m 4 Si 4 S24 - ((mi - m 2 )s 12 c 12 c 23 - m 3 s 23 s 13 ^J = , (57d) 

m er : m A s u s 3A + ((mi - m 2 )s 12 c 12 s 23 + m 3 c 23 Si 3 ) = , (57e) 

m MT : m 4 s 2 4S 3 4 - S2 3 c 23 (mis 42 + m 2 c\ 2 - m 3 ) = . (57f) 



Obviously, the first four equations can never be satisfied, as the left-hand sides of them 
are definitely positive. Consequently, none of the 20 type II textures without CP violation 
can even be compatible with the current data on active neutrino oscillations. We calculate 
numerically all these 20 CP conserving textures and predictively obtain the negative result. 
For the 100 type I CP conserving textures, due to the same reason above for the type II tex- 
tures, they are all expected to be excluded by the active neutrino data. We also numerically 
analyze these 100 CP conserving textures and again arrive at the same conclusion. 



VI. CONCLUSION 

The LSND anomaly, MiniBooNE experiment, the reactor anomaly, the Gallium anomaly 
and cosmological observations all point to the existence of light eV scale sterile neutrino(s). 
Texture zero is a powerful framework to investigate the phenomenologies of the neutrino 
sector of the SM. In this work we apply the framework of texture zeros to the smallest 
extension of the neutrino sector of SM, i.e., the (3+1) model, and examine all the 120 
textures of 4 x 4 Majorana neutrino mass matrix M v with three zeros. Although all the 
3x3 neutrino matrices with three zeros have been excluded, we find that when the neutrino 
sector is extended to include one sterile flavor, the situation is completely different. 

We classify the 120 textures into two large classes, the 100 type I textures with zeros in 
the sterile sector (the fourth row and column) of M v and the 20 type II textures without 
zero in the sterile sector. With the approximations of large mass hierarchy 2,3^ ^4 and 
smallness of mixings #13,14,24,34 *C 1, all the analytical expressions can be greatly simplified. 
We analyze both analytically and numerically all the 120 textures and find that the 100 type 
I textures are all excluded by current neutrino data, mainly due to the large active-sterile 
mass splitting. The exception among the type I textures are the class £ textures with only 
m TS = in the sterile sector, but they are highly disfavored by experiments as a result of 
the extremely small available parameter space. The 20 type II textures with all the three 
zeros in the active sector are phenomenologically much more favored. As for the textures 
within the three- neutrino scheme [231 ]. some of the textures are related by the approximate 



fi — r symmetry [25[ and the symmetric partners always share the same phenomenological 
features. Their remarkable features revealed by both the analytical and numerical analysis 
are collected as follows: 

1. To obtain vanishing elements in the active sector of M v and keep compatible with the 
(9(0.1) active-sterile mixing, the active neutrino mass and lightest neutrino mass are 
required to be of order 0.01 eV. 

2. We can also state the same fact in another way: If the active neutrinos are of order 
0.01 eV, we naturally arrive at the ASMs are of order 0.1. That is to say, the 20 
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textures are overall compatible with current experimental data, which is verified by 
both analytical and numerical analysis. Numerical calculations reveal further that 
among the 20 textures, only the texture C3 is actually excluded and the other 19 are 
more or less compatible with data. 

3. Both the two neutrino spectrums of NH and IH for the three light active flavors are 
allowed in the textures A, 62,4,5,6, D 3 4 5 6 , E 2 ,3, F, whereas the seven textures 
Bi,3, C12, Di,2, Ex prefer much for the NH spectrum over IH. 

4. Some of the CP violating phases in the 19 textures exhibit preference of particular 
values, e.g., the Dirac phase 5i 3 in the texture F favors the values around n/2 or 3n/2, 
and in A the Majorana phase a ~ n. All the predictions on the phases are collected 
in Table [VTl The Dirac- type phases, especially the phase #13, can be observed in the 
upcoming long baseline neutrino oscillation experiments, while the Majorana phases 
are crucial to the OvflP experiments. All these relevant oscillation and nonoscillation 
experiments are promising to examine the variant phenomenological implications and 
test these textures in the (3+1) model. 

5. The 0.01 eV active neutrinos and eV scale sterile neutrino sum up to ~eV, compatible 
with the cosmological bound on neutrino mass. No matter whether the KATRIN beta 
decay experiment can observe the electron neutrino mass or not, it is challenging for it 
to distinguish the two neutrino spectrums in the three- zero textures. In contrast, it is 
promising for the next generation 0^/3/3 experiments with sensitivity of order 10 meV 
can clearly differentiate NH from IH and then test these textures. For example, if the 
effective Majorana neutrino mass comes up to be 30 meV, the three active neutrino 
are probably of inverted hierarchy and then those textures preferable of NH would be 
almost excluded. 5 

In addition to the textures with the most general CP violation, we calculate also the textures 
without any CP violation and find that without the tuning of the phases in the correlation 
equations, none of the 120 textures are allowed by current neutrino data. 

In the end we stress that with more accurate neutrino data, especially the determination 
of the octant of #23 1381, the textures with three zeros in the (3+1) model will be more 



severely constrained. It is likely that more textures will be excluded and we can obtain more 
accurate phenomenological predictions, e.g., the effective neutrino mass in the neutrinoless 
double beta decay. 
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5 The 0v(3f3 test of the textures in this item is based on the effective neutrino mass, while the Oi//3/3 test in 
the previous item is related to the Majorana phase detection in Oi//3/3 experiments. Although both the 



tests are in principle feasible, the latter one is much more challenging than the former 
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Appendix A: Classification of the 120 textures 



The number of textures for 4x4 M v with three zeros is rather large; thus, it is beneficial 
for us to properly classify and label the 120 distinct textures. As the sterile neutrino is 
essentially different from the three active ones, the 120 textures can be roughly classified 
into two large categories: (I) 100 textures with zeros in the sterile sector of M v and (II) 20 
textures without vanishing elements in the sterile sector. As we will see in the analytical and 
numerical discussions, predictions of the structures in the two classes on the ASM angles 
are completely different. 

For the convenience of discussions, the 100 textures of M v of class I can be further 
classified into six subgroups A — J 7 , according to the number of zeros in the sterile sector. 
The matrices given below are the sample matrices for the six subgroups: 
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All the three texture zeros of the first two groups A and B are in the sterile sector. The 
three textures of type A have the fourth diagonal element m ss = 0, whereas all the three 
ASM angles of the texture of type B are zeros and this texture is apparently excluded by 
the LSND experiment. The two following groups C and T> consist of the M v structures with 
one zero in the active sector and the other two in the sterile sector. m ss = holds true for 
the textures of type C, while for the textures of type T>, the two texture zeros in the sterile 
sector are both nondiagonal elements. For the last two groups £ and J 7 , only one of the 
three texture zeros are from the sterile sector. For type £ m ss = 0, whereas for type T one 
of the three nondiagonal elements in the fourth row or column vanishes. 

The 20 patterns of M v with all the three zeros in the active sector can also be further 
classified into subgroups. 3x3 neutrino matrices with three texture zeros have been classified 
according to the pos itions of zeros, or more specifically according to the number of vanishing 
diagonal entries [18J. The 20 4 x 4 matrices are analogically classified into six subgroups 
and they are explicitly labeled as follows: 
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/ x x x \ 

x x x 

x x x 

x x x x 



(A2) 



/ x x \ 

x x x 

x x x 

x x x x 



/ x x x \ 

x x 

x x x 

X X X X 



B. 
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(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



All the first three diagonal elements of the texture A vanish, whereas two of the first three 
diagonal elements of the textures of types B and C are zeros, only one for types D and E 
are vanishing and none of the three diagonal elements in the active part of F are zero. The 
three textures of type C are of rank three, which are distinguishable from the six of type B 
which are of rank four. Similarly, the six structures of type D are of rank four, whereas the 
three of type E are of rank three. 
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